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We investigate the combined effects of magnetic impurities and applied magnetic field on the 
interference contribution to the conductance of disordered metals. We show that in a metal with 
weak spin-orbit interaction, the polarization of impurity spins reduces the rate of electron phase 
relaxation, thus enhancing the weak localization correction to conductivity. Magnetic field also 
suppresses thermal fluctuations of magnetic impurities, leading to a recovery of the conductance 
fluctuations. This recovery occurs regardless the strength of the spin-orbit interaction. We calculate 
the magnitudes of the weak localization correction and of the mesoscopic conductance fluctuations 
at an arbitrary level of the spin polarization induced by a magnetic field. Our analytical results for 
the "ft/e" Aharonov-Bohm conductance oscillations in metal rings can be used to extract spin and 
gyromagnetic factor of magnetic impurities from existing experimental data. 



I. INTRODUCTION 

Conductance of disordered metals is sensitive to the 
applied magnetic field. At weak magnetic field the con- 
ductivity of a bulk metal has a sharp peak due to the 
weak localization (WL). 1 Similarly, the conductance of 
mesoscopic metals fluctuates as the magnetic field is 
changing. 2,3 Both the conductance fluctuations and the 
WL correction to the conductivity are quantum mechani- 
cal phenomena originating from the interference of quan- 
tum states. As any other interference phenomena, they 
may be demolished by interaction processes. 4 

Localized spins affect the electron transport in metals. 
Various properties of electron kinetics are sensitive to dif- 
ferent aspects of the localized spin dynamics. The energy 
exchange between electrons in the process of scattering 
off a magnetic impurity is made possible by the quan- 
tum fluctuations of the impurity spin: its virtual flip in 
the course of scattering facilitates the energy transfer be- 
tween the two electrons. 5 On the other hand, no spin 
dynamics of impurities is needed for the suppression of 
the weak localization correction to conductivity; interac- 
tion of electron spins with randomly oriented magnetic 
moments already leads to that suppression. 6 

Universal conductance fluctuations (UCF) are not sup- 
pressed by static magnetic moments. However, even 
a relatively slow relaxation of individual magnetic mo- 
ments, such as Korringa relaxation, 7 leads to the time- 
averaging of the random potential "seen" by transport 
electrons in the course of measurement, and the meso- 
scopic fluctuations of the dc conductance get averaged 
out. 8,9 The sensitivity of conductance fluctuations to the 
time evolution of the system of localized magnetic mo- 
ments 10 was used extensively to probe the spin-glass 
freezing in metals 11,12 and semiconductors. 13 

An applied magnetic field may be used to control the 
statistical properties and dynamics of localized magnetic 
moments. It was noticed 12 that in a strong magnetic field 
the amplitude of conductance fluctuations and "h/e" 
Aharonov-Bohm (AB) oscillations increases, apparently 
because the spin fluctuations are quenched. 8 Recently the 
dependence of the amplitude of "/i/e" AB oscillations on 



magnetic field was carefully measured on Cu wire rings. 1 
The goal of that measurement was to corroborate the 
existence of localized spins, conjectured on the basis of 
measurements of the electron energy relaxation rate in 
Cu wires. 15 

Theory of conductance fluctuations and WL correction 
to the conductivity at partial spin polarization of mag- 
netic impurities has not been developed yet. Only the 
limits of no spin polarization at B = and of strong po- 
larization at B ^> T/gfiBS were considered 8, 16-18 (here 
B is the magnetic field, S is the impurity spin, T is the 
system temperature, g is the impurity gyromagnetic fac- 
tor, and fj,B is the Bohr magneton). 

In this paper we concentrate on the interference contri- 
bution to the linear conductance of mesoscopic systems in 
the presence of partially polarized magnetic impurities. 
Particularly, we calculate the weak localization correction 
to the conductivity and the amplitude of the mesoscopic 
conductance fluctuations. 



II. MAIN RESULTS 

It is well known that scattering of electrons off mag- 
netic impurities in the absence of a magnetic field sup- 
presses the interference correction to the conductivity of 
a wire: 1,6 

Aa = __ !!_VTJ ( -JL= - — = V (l) 

Here D is the diffusion constant for electrons in the wire 
and 1 /r s is the electron scattering rate off magnetic im- 
purities. To the lowest order in the exchange constant 
J calculation, l/r s = 2nvn s J 2 S(S + 1), where n s is the 
concentration of magnetic impurities, v is the electron 
density of states at the Fermi level per spin degree of 
freedom, and S is the magnitude of impurity spins. Tak- 
ing into account the Kondo renormalization of the ex- 
change constant, at temperatures well above the Kondo 
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temperature Tk, this rate can be cast in the form 
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It is clear from Eq. (2) that the phase relaxation rate 
increases upon reduction of the temperature towards Tk- 
This rate reaches maximum at T <~ Tk, and decreases if 
the temperature is further reduced. The specific form 
of function r s (T) at low temperatures depends on the 
spin 5 of local moments. The screening of the local spin 
5 = 1/2 is complete at T — and at T <C Tk the phase 
relaxation rate can be found from the Nozieres' Fermi 
liquid theory, 
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If S > 1, then screening is incomplete, and the loga- 
rithmic renormalization theory can be employed for the 
evaluation of r s at low temperatures, 19 
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Note, that the low-temperature asymptotes, Eqs. (3), dif- 
fer from the approximate expression, used sometimes in 
the analysis of experiments. 20,21 

In this section, we only present the result for WL cor- 
rection Act to the conductivity of a thin wire in the ab- 
sence of spin-orbit scattering. If, in addition, the density 
of magnetic impurities is small, then starting from rather 
weak fields the condition g^ B Br 8 » 1 is satisfied, and 
we find 
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Here A is the wire cross-section area, i? is a numerical 
factor, and <&o = hc/2e. Function F(e) represents the 
spin flip rate in the presence of magnetic field B: 



T(e) 
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Here (...) stands for the thermodynamic average over 
the states of an isolated impurity spin, see Eq. (15). 
The Kondo-renormalized rate 1/t s in Eq. (5) is given by 
Eq. (2) at weak magnetic field, B < T/gfiB, and by 
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at strong magnetic field, B > max{T, T^/g^B- 

In the domain of relatively weak magnetic fields, 
1/g^BTs <C B <C T/g^B, the spins arc not polarized, 
however, their precession already affects the value of WL 
correction. Equation (4), valid in this domain, yields 
Act differing from the full zero-field value (1) by a factor 



2/(3 — \/l/3) ~ 0.83. We discuss the detailed behavior 
of the WL correction in the field region B < l/gf^BT s 
in Section IV. The main variation of the WL correc- 
tion, however, is associated with the spin polarization. 
This variation occurs at B <~ T/ g^B , and is described by 
Eqs. (4) and (5). 

At strong magnetic field, B > T/gfiB, the spin flip 
rate of thermal electrons (e <~ T) becomes exponentially 
small, r cx cxp(— g^B /T) /[(S + 1)tJ. In this case 
higher-order terms in the exchange interaction must be 
retained to calculate the phase relaxation rate. We find 
the following result for T(e): 



T(e) 
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The crossover between the rates given by Eqs. (5) and 
(7) occurs at 
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The relaxation rate at high fields, Eq. (7), may still 
be large enough to compete with the conventional 
mechanisms 1,4 of phase relaxation caused by the electron- 
electron and electron-phonon interactions. 

If one disregards the orbital effect of magnetic field 
(very thin wire) then the WL correction to the conduc- 
tivity, Eq. (4), is a function of the ratio B/T, up to slowly 
varying logarithmic factors \n{gix B B/TK) and Iii(T/Tk), 
see Eqs. (2) and (6). 

Different regimes for the WL are schematically rep- 
resented in Fig. 1. Substantial variation of Act with B 
occurs at B > T/g^S- The curve B = B*(T) is shown 
in Fig. 1 by the bold solid line. 

The effect of magnetic-field-induced polarization of lo- 
calized spins on the WL correction to the conductivity is 
significant only if two quite stringent conditions are met. 

First, the orbital effect of the magnetic field B on WL 8 
is needed to be small in the range of fields B > T/g/iB, 
where the spin polarization is significant. The orbital 
effect dominates over the spin scattering, if DAB 2 /<&q < 
r(e ~ T), see Eq. (4). This condition defines the upper 
bound (a dashed line in Fig. 1) for the domain of fields B 
in which the spin polarization noticeably affects the WL 
correction. In general, to avoid the undesirable orbital 
effect, measurements are to be done on thin wires at low 
temperature. 

Second, the spin-orbit scattering must be sufficiently 
weak. This condition limits the range of host conductors 
to materials with sufficiently low atomic numbers. 1 In 
heavier materials, the spin-orbit interaction suppresses 
the triplet part of the WL correction to the conductiv- 
ity, represented by the first term in Eq. (1). Polarization 
of local moments does not eliminate the randomness of 
the effective magnetic field induced by spin-orbit interac- 
tion, and the triplet contribution to the WL remains sup- 
pressed. At the same time, the singlet part of the weak 
(anti)localization correction [the second term in Eq. (1)] 
is insensitive to the spin polarization. 
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FIG. 1: Different regimes of the temperature, T, and mag- 
netic field, B, dependence of the mesoscopic conductance 
fluctuations and the WL correction to conductivity. Below 
the Kondo temperature, T<Tk, and at weak magnetic field, 
qUbB^Tk , magnetic impurities are screened and the elec- 
tron phase relaxation rate is given by Eqs. (3). The impurity 
magnetic moments acquire a significant polarization at fields 
around B = T/gfxsS (straight solid line), which results in 
strong dependence of the electron phase relaxation rate on 
the applied magnetic field. The exponential behavior of the 
phase relaxation rate on B/T is replaced by a weaker power- 
law function, see Eq. (7), at a crossover field B* defined by 
Eq. (8) and depicted by the solid bold curve. Strong varia- 
tion of the WL correction is possible only in the absence of 
spin-orbit interaction and at sufficiently weak orbital effect of 
the magnetic field. This last condition is satisfied below the 
dashed curve B a (T) defined by Eq. (34), which intersects the 
line gusSB = T at temperature T Q . 



The effect of spin polarization on the mesoscopic con- 
ductance fluctuations is not subject to the two restric- 
tions discussed above. In other words, the manifestation 
of the spin polarization in the conductance fluctuations 
and in the u h/e" Aharonov-Bohm effect is much more 
robust than that in the WL. The presence of magnetic 
impurities suppresses the conductance fluctuations only- 
due to the spin dynamics of the impurity system. If the 
measurement time is significantly longer than the char- 
acteristic time of the variation of impurity spin config- 
uration, then the conductance fluctuations are averaged 
out. Quenching of the spin configuration by the applied 
magnetic field decreases the effect of the conductance 
averaging and restores the fluctuations. In the limit of 
full polarization, magnetic impurities no longer affect the 
conductance fluctuations. 8 

We concentrate on the amplitude of the Aharonov- 
Bohm "/ic/e" conductance oscillations, since they are 
exponentially sensitive to the polarization of impurity 
spins. Magnetic field flux $ threading the ring of radius 



R changes electron wave functions and, consequently, the 
conductance g$> of the ring. The conductance statistics 
is characterized by the correlation function: 
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where $ = irR 2 B is the magnetic flux through the ring. 

We find the amplitude of oscillations of the conduc- 
tance correlation function: 
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where 



T(e) = 



1 - 



(S z ) 2 + (S z ) tanh(e + g^ B B)/2T 



S(S + 1) 



(11) 

is the phase relaxation rate to the second order in the 
exchange interaction J, and a is a dimensionless geom- 
etry dependent factor. As the magnetic field increases, 
r(e) decreases exponentially and at field given by 
Eq. (8), the higher order term in the exchange interac- 
tion, Eq. (7), prevails. 

Equation (10) is valid at T(e ~ T) « T. This condi- 
tion is always satisfied at a sufficiently strong magnetic 
field, since T(e ~ T) — > as B increases and the impu- 
rity spins become almost polarized, see Eqs. (7) and (11). 
Here again we disregarded other mechanisms of phase re- 
laxation, such as the electron-electron interaction. 1 ' 22 

We emphasize that Eq. (11) contains the average (S z ) 2 
instead of the term (S 2 ) in Eq. (5). The two equations 
are different because of the nature of the electron dynam- 
ics producing the WL correction to the conductivity and 
conductance fluctuations. The difference can be under- 
stood in the following way. 

The weak localization correction to the conductiv- 
ity originates from the electron passage along the same 
trajectory twice. The time difference between the two 
passages does not exceed the time of phase relaxation, 
[r(e ~ T)]~ x in our case. In the derivation of Eq. (5) we 
assumed this time being shorter than the Korringa re- 
laxation time 7 tt, see Eq. (18). Hence the instantaneous 
spin-spin correlator {S 2 ) enters Eq. (5). The origin of 
the correlations in mesoscopic conductance is also due 
to passage of electrons along the same trajectory. Here, 
however, the relevant time difference is defined by the 
time between the measurements, and typically greatly 
exceeds tt- That is why the conductance correlation 
function is described by (S z ) 2 correlator, characterizing 
a non- fluctuating component of impurity spins. 

Note that in the case of low concentration of magnetic 
impurities n s < vT the Korringa relaxation time may 
become shorter than [T(e <~ T 1 )] -1 . In this case Eq. (11) 
rather than Eq. (5) defines the phase relaxation rate in 
the WL correction to the conductivity. The correspond- 
ing modification of the WL correction at B = was 
considered earlier in Ref. [18]. 
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To summarize, wc studied the effect of magnetic field 
on the weak localization correction and on the magnitude 
of mesoscopic conductance fluctuations in a conductor 
with magnetic impurities. Our results are valid at an 
arbitrary level of the spin polarization. We demonstrate 
that the electron phase relaxation rate acquires energy 
dependence due to the Zeeman splitting. Such energy de- 
pendence is absent in the limits of no impurity polariza- 
tion, B = 0, and strong polarization, g/isB 3> T. How- 
ever, in the experimentally relevant intermediate regime, 
gfJ.sB ~ T, the effect of spin polarization can not be ac- 
counted for by assigning a single phase relaxation rate to 
all electron states. 



III. MODEL 

We consider a metal with isotropic elastic scattering of 
electrons by impurities, characterized by the mean scat- 
tering rate of electrons l/r c = . Here v is the 
electron density of states at the Fermi surface per spin 
state, n c is the concentration of impurities, and U is the 
Born amplitude of elastic scattering by an impurity. 

The scattering of conducting electrons off magnetic im- 
purities is described by the Hamiltonian: 



H m = JSa, 



(12) 



where 5 is the spin operator of a magnetic impurity, and 
J is the exchange constant. The electron scattering rate 
by magnetic impurities is 
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Here n s is the magnetic impurity concentration, and 
5 is the total impurity spin. The exchange constant 
J is renormalized due to the Kondo effect: J — > 
2/ '{v\xiT /Tk) at temperatures T exceeding both Tk and 
gH&B. In stronger fields, B > max{T, Tk}/ g^B, tem- 
perature T under the logarithm is replaced by <?/kb-B. 

We study the effect of magnetic impurities on conduc- 
tivity of metals. In order to evaluate the interference 
effects in electron transport we will need instantaneous, 
(S z ) and (5 2 ), and time delayed spin-spin correlation 
functions: 



Xz(t) 
X±(r) 



(S z (t)S z (t - t)) 
5(5+1) ' 
([5 + (t),5_(f-r)] + ) 
5(5 + 1) 
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Here S z (t) and S±(t) = S x (t) ± iS y (t) are parallel and 
perpendicular spin components of a magnetic impurity, 
[■,■]+ is an anticommutator, and 
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stands for the thermodynamic average in the presence 
of a magnetic field B, producing the Zeeman splitting 
w s = g^B,B of the magnetic impurity states; g is the 
impurity spin gyromagnctic factor. 

For a metal with dilute non-interacting magnetic im- 
purities, functions Xz( T ) and XJ-( T ) have the form: 

, , W + {(5 2 ) (5 2 ) 2 } f 2 (r) 
Xz(t) = g . p ■ 7T , (16a) 



X-l(t) = 



5(5 + 1) 
{5(5+1)-(5 2 )}/ ± (t) 
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The functions f z (r) and f± (t) describe spin relaxation. 
In a dilute system of magnetic impurities, relaxation oc- 
curs due to the Korringa mechanism 7 and is exponential: 



f±(r) 
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Here tj_ and t z are the spin relaxation times for compo- 
nents perpendicular and parallel to the applied magnetic 
field, respectively. In a weak magnetic field the two spin 
relaxation times coincide and are equal to the Korringa 
time tt- 
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The applied magnetic field produces Zeeman splitting 
£z of energy states of transport electrons. When esti- 
mating the Zeeman energy ez in metals with magnetic 
impurities, it is important to notice that the sign of the 
exchange constant J is fixed by the nature of the pair of 
host and impurity 13 atoms. Polarization of magnetic im- 
purities results in the exchange contribution n s J(S z ) / fiB 
to the effective magnetic field which causes the Zeeman 
splitting £z of spin states of transport electrons with the 
gyromagnetic factor g e 



ez = g c m B - 2n s J(5 z ). 



(19) 



Equation (19) shows that magnetic impurities may signif- 
icantly affect the Zeeman splitting of transport electron 
states. 

In the case of antiferromagnctic exchange, J > 0, the 
dependence of £z on B is not monotonic at low temper- 
atures. As we will see below, Eq. (19) provides a mecha- 
nism for a non-monotonic in B interference contribution 
to the metal conductivity, similar to the Jaccarino-Peter 
mechanism of the reentrant superconductivity. 23 ' 24 

We also consider the effect of spin-orbit scattering on 
the conduction of electric current. For the specific case 
of heavy element impurities of concentration n so , the 
corresponding term of the Hamiltonian has the form: 25 
H so = U so [p x p']a, and the spin-orbit scattering rate 

7so = 27wn S0 L r2 ?4- 

In this paper we assume that the shortest electron re- 
laxation time is due to elastic scattering and calculate the 
conductivity to the lowest order in t c "/ s and r 7 so , using 
the standard diagrammatic technique for a disordered 
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metal. 26 Nevertheless the results are valid even for con- 
ductors with strong spin-orbit interaction in their host 
material, when the spin-orbit scattering rate is compara- 
ble with elastic scattering rate. 



IV. WEAK LOCALIZATION CORRECTION TO 
THE CONDUCTIVITY 

The weak localization correction to the conductivity is 
given by the following expression: 



e 2 D 
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where the Cooperon C(e,w, q) is 



(20) 



(21) 



The Cooperon components, corresponding to different in- 
dices i, represent possible spin configurations: indices 
i = 1,2 are assigned to the Cooperon components with 
non-zero spin projections on the direction of magnetic 
field, m c = ±1 and S c = 1; index i = 3 is assigned 
to the Cooperon with S c = 1 and zero spin projection 
m c = 0; and index i = 4 stands for a singlet component, 
S c = m c = 0, sec Table I. 

Performing integration over momentum q, we obtain 
a general result in the form 



(22) 



Functions ^(B, T) have different form depending on the 
conductor geometry (d — 2 for a metal film and d = 1 for 
a wire), while parameters oti and phase relaxation rates 
Fi(B) are determined by details of electron scattering 
processes. First we present functions JFd(B,r) for differ- 
ent conductor geometries and then discuss the effect of 
spin flip scattering on the weak localization correction to 
the conductivity. 



A. Geometry dependence 

The weak localization correction to the conductivity 
of a quasi two dimensional metal is determined by the 
following expression: 



Fi (B,Tj) 
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47T 



A I\ + 7 orb 
P \2 ADeB ± /hc 



+ In 



AeB ± Dt c 

he 



(23) 



TABLE I: The Cooperon phase relaxation factors for various 
components corresponding different spin configurations. The 
components with i = 1, 2, 3 represent a spin states with total 
spin one, Sc = 1- These spin configurations form a triplet in 
an isotropic system. Anisotropy due to the applied magnetic 
field splits the degeneracy. The i — 4 component is a singlet 
spin configuration with zero total momentum, Sc = 0. 



i 


!Sc,m c ) 






1 

2 
3 
4 


S c = l,m c = +1 
S c = l,m c = -1 
S c = l,m c = 
Sc = 0, mc = 


i - (Si)/s- 
i - (sl)/s 2 

2<S 2 }/S 2 
2 


+1 
+1 
+1 

-1 



where ip(x) is the digamma function and I\ is the dephas- 
ing rate for the ith component of the Cooperon. Here Tj 
includes dephasing 7 or b by the applied magnetic field par- 
allel to the film 7 or b = e 2 £?jj Da 2 /12h 2 c 2 , where Bj_ and 
By are the magnetic field components, perpendicular and 
parallel to the film. In weak perpendicular magnetic field 
we use '0(1/2 + x) s=a h\x + l/24x 2 and obtain: 



^(B.ro 



i 

47T 



In- 



1 



(r, + 7orb)r c 



2(DeBj_/hc 



3(r.t + 7orb) 2 



(24) 



The terms with T e in Eqs. (23) and (24) originate from 
the ultraviolet logarithmic divergence of the momentum 
integral for the weak localization correction to the con- 
ductivity in two dimensional case. The existence of this 
term complicates analysis of the magnetoresistance mea- 
surements. 

In the quasi one dimensional case the integral over mo- 
mentum gives: 



^i(B,r 4 ) = 



D 



1 



2 V + 7orb 



(25) 



The denominator of Eq. (25) contains the sum of electron 
phase relaxation rates due to scattering on magnetic im- 
purities and due to the orbital effect of the applied mag- 
netic field 7 or b- For a magnetic field of strength B and 
applied in direction n, the orbital contribution to the 
phase relaxation rate is 



7 orb = i?(n) 



e 2 B 2 DA 
h 2 c 2 ' 



(26) 



where the function i?(n) of the direction n is of the order 
of unity, and A is the area of the wire cross-section. 

In the following subsections we focus on properties of 
the phase relaxation rate r, due to scattering on mag- 
netic impurities. To be specific, we will present the WL 
correction to the conductivity of a thin wire. Since the 
phase relaxation factors 1^ do not depend on the sam- 
ple geometry, the results for the WL correction to the 
conductivity of the next two subsections can be easily 
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FIG. 2: Dependence of the weak localization correction on 
the applied magnetic field in a metal with classical spin im- 
purities. The plot shows function Y(B, T) at two temperature 
values T = 0.3T o and T = 0.03T o for both short (r T < r s ) and 
long (tk 2?> t s ) impurity spin relaxation time tt- The sub- 
plot shows the WL correction to the conductivity for S — 1/2 
impurities at T — T . Here T is given by Eq. (35). 



and for the phase relaxation rates 

Ti,2 = (1-X*)-, ( 2 8a) 
T 3 = {l + Xz)-+ius-Z, (28b) 
T 4 = {l + x*)-+ius + Z. (28c) 

T s 

Here 




- (e z +ws) 2 ; 



Xz and xj_ are the impurity spin correlators given by 
Eqs. (16) and (17) with proper substitution r = or r — > 
oo. The correlation functions (S z ) and {S%} of classical 
spins are described by the following functions: 

(S z ) = Scoth^-— , (29a) 



generalized to other system geometries, such as a film, a 
metal ring or an open quantum dot. 27 



The weak localization correction to the conductivity of 
a wire may be represented as 



B. Effect of classical spin impurities 

Even in the absence of spin-orbit interaction, there are 
five parameters affecting the WL correction to the con- 
ductivity, including system temperature T, the applied 
magnetic field B, impurity spin relaxation time tt, the 
electron scattering rate on magnetic impurities 1/t s , and 
the phase relaxation rate due to the orbital effect of the 
applied magnetic field 7 or b- It is easier to establish the 
role of these parameters in the case of classical (S> 1) 
spins, considered in this Subsection. 

The explicit matrix form of the equation for the 
Cooperon in the case of classical spins is given in Ap- 
pendix A, see Eqs. (A2) therein. For an arbitrary relation 
between the electron dcphasing rate Ti and the impurity 
spin relaxation rate 1/tt, the solution of this equation 
is cumbersome. Here we consider only the limiting cases 
of long (r^Ti >• 1) and short (tiT, <C 1) spin relaxation 
time. Usually (for not too small concentration of mag- 
netic impurities n s > vT) the first condition is satisfied; 
the second condition may become relevant only at very 
low n s . In the first case one may neglect 6 the dynamics of 
localized spins setting t = in the correlation functions 
Eq. (A2b). In the second case, we need to set r — > oo 
in Eq. (A2b). In both limits we obtain the following ex- 
pressions for the parameters of Eq. (22) 



ai 2 = 1, a.3 = + 



X-L 

Zt s 



Q!4 



X-L 



(27) 



Aa 1 = - — VDr s Y(B,T). 



(30) 



The function Y(B,T) has different forms, depending on 
the relations between parameters of the system. 

In the limit of long spin relaxation time (rrTi 3> 1) 
and weak magnetic field (wst s <C 1 and £zt s -C 1) all 
four Cooperon modes contribute to the weak localization 
correction to the conductivity 



Y(B,T)= 



+ 



v/l-(,S 2 )/5 2 +7orbTs 

1 1 



y/2(S*)/S 2 +'y olb T a V2 + 7orbT 5 



• (31) 



The first term in Eq. (31) is due to the Cooperon modes 
with i = 1,2. It is the first term of the function Y(B,T) 
which is responsible for a non-monotonic magnetic field 
dependence of the WL correction to the conductivity. 
As magnetic field increases and impurity spins become 
polarized, the phase relaxation rate for modes i = 1,2 
decreases, and the first term in Eq. (31) grows. The 
combination of the remaining two terms, representing 
Cooperon components i = 3,4, constitutes 0.17Acri of 
the full value of Acti at zero magnetic field and mono- 
tonically decreases as the magnetic field increases. 

We notice that the contribution of i = 3,4 Cooperon 
components diminishes not only as magnetic impurities 
become polarized, but also as the Zeeman splitting of 
conduction electrons and impurity spins (ez + ws) m- 
creases, see Eqs. (20) and (27). At strong magnetic 
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field, \u>s + £z\i~ s S> 1, the contribution of i = 3,4 
Coopcron components to Aai vanishes, since both the 
Zeeman splitting of the conduction electrons and the en- 
ergy transferred in the spin flip process destroy the in- 
terference of time reversed paths in the Cooperon. If the 
spin relaxation time is still long, "rrTj » 1, we have 



Y (B, T) 



v/l-(S2)/S 2 + 7 orbT B 



(32) 



For impurity spin system with short spin relaxation 
time, Trl\ <C 1, the Cooperon modes with i = 3, 4 do not 
contribute to Acti again. In this case, the weak localiza- 
tion correction to the conductivity is given by Eq. (30) 
with 



Y(B,T) = 



V /l-(,S z ) 2 /5 2 +7orbT s ' 



(33) 



We notice that as magnetic field increases, the phase 
relaxation rate decreases, and the condition of short 
spin relaxation time (ttTi^ <C 1) may be reached, 
even though initially system was in the opposite regime, 
tt > r a . 

We conclude from Eqs. (31)-(33) that the WL cor- 
rection to the conductivity depends on the polarization 
of magnetic impurities through the thermodynamic av- 
erages (S z ) and {Si). The polarization is significant in 
an applied field of the order of B = T/gfi^S, shown as a 
straight solid line in Fig. 1. 

The orbital effect of the applied magnetic field is repre- 
sented by the term 7orb7s i n Eqs. (31)-(33). The orbital 
effect and the impurity spin polarization compete with 
each other: the orbital effect suppresses, while the spin 
polarization enhances the WL correction to the conduc- 
tivity. The phase relaxation rates 7 or b and Ti become 
equal, 7 or b = T\, at magnetic field B (T): 



Bo(T) 



VDt~A 



1 - 



(S z )%-_ 



B 



s 2 



(34) 



The effect of spin scattering prevails over orbital effect 
of the magnetic field at B < B Q (T). On the other hand, 
the enhancement of the WL correction due to impurity 
spin polarization is significant at B > T/g^S. These 
two limitations are met if the wire temperature T < T , 
where 



T = g^ B S 



(35) 



At higher temperature, T > T , the effect of the spin 
polarization is concealed by the orbital effect. 

In Fig. 2 we illustrate the effect of an applied magnetic 
field on the WL correction to the conductivity. We show 
dependence of the function Y(B, T) on the applied mag- 
netic field B at temperatures T = 0.3T o and T = 0.03T o . 
At low temperature T — 0.03T o , the effect of impurity 
spin polarization is significant, but already at tempera- 
ture T = 0.3T o it fades away. 



The spin-orbit interaction modifies the WL correction 
to the conductivity. Using Eqs. (A2) in Appendix A we 
find: 



Oil, 2 = 1, 



a 3 



27so/3 - X-L/f s 



V(2 7 so/3-x±M) 2 -^ 



Q!4 = — - 



27so/3 - X_l/t s 



v /(2 7 so/3 - X ±/r s ) 2 - e 
and the phase relaxation rates are 

ri,2 = (l-Xzz) 

1 + Xzz , 27 



(36a) 
(36b) 

(36c) 



(37a) 



r, = 




-si (37b) 



l + Xzz, 27s, 
1 4 = 1 — — 



X± 2 7s, 
t s 3 



e\. (37c) 



Here we assumed us <C 1^ for simplicity. 

According to Eqs. (37), in metals with strong spin- 
orbit scattering, 7 so t s 3> 1, the terms with i — 1,2,3 
are suppressed and the interference correction to the con- 
ductivity is described by the singlet antilocalization term 
i — 4. In the limit of long impurity spin relaxation time 
tt ^S> t s , the weak localization correction to the conduc- 
tivity depends on r s and the Zeeman energy ez- 



Y(B,T) = 



1 



y/2 + 3e|r s /47 so + 7 orb r s ' 



(38) 



while in the opposite limit tt <§C t s the correction also 
depends on the polarization of impurity spins: 



Y(B,T) = 



y/1 + (S z ) 2 /S 2 + 3£ 2 r s /4 7so + 7orb r s 



.(39) 



The terms 7 or b7"s oc A in the denominator of Eqs. (38) 
and (39) originate from the orbital magnetic field effect, 
see Eq. (26). If the area A is sufficiently small so that 
the orbital part is not important, and the temperature is 
low enough, then the antilocalization may show a non- 
monotonic field dependence due to the Jaccarino-Peter 
mechanism, 23 see Eq. (19). 

From the above analysis performed in the approxima- 
tion of large impurity spins we make two conclusions. 
First, the polarization of impurity spins hardly affects 
the weak (anti)localization correction to the conductiv- 
ity in metals with strong spin-orbit interaction. Indeed, 
the only surviving contribution to the conductivity cor- 
rection originates from the singlet spin configuration of 
the Coopcron, which is not sensitive to the impurity spin 
polarization. Second, in metals without spin-orbit inter- 
action, the Cooperon modes with zero spin projection on 
the magnetic field (i = 3, 4) give only a small (~ 0.17Acr) 
contribution to the total WL correction to the conductiv- 
ity, and their contribution vanishes at moderately strong 
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FIG. 3: The Cooperon decay rate is determined by diagrams of the second order in the exchange constant J. The two left 
diagrams represent the electron self energy contribution from processes with and without (not shown) spin flip. The right 
diagram represents the vertex correction. 



(uJsT s ~ 1) magnetic field. In the next subsection we ne- 
glect these terms and calculate the WL correction to the 
conductivity, originating from Cooperon modes i = 1 , 2 
at such fields that ust s > 1. We perform the calculations 
for an arbitrary value of S. 

C. Effect of quantum spin impurities 

We notice that the semiclassical description is applica- 
ble only i) for a large spin S » 1, when linear in S con- 
tributions to the dephasing rates can be neglected and 
ii) at high temperature T 3> us, so that the discreetness 
of spin energy levels can be disregarded. In this Sub- 
section we consider quantum spins with arbitrary value 
of S and for arbitrary ratio of us/T. As we discussed 
in the previous Subsection, only components i = 1,2 of 
the Cooperon with parallel spins, see Table I, are impor- 
tant at strong magnetic field. The terms with i = 3, 4 are 
small at relatively strong magnetic fields, \us+£z\t s > 1. 
Moreover, even at weak fields, modes i = 3, 4 give rise 
only to a small portion of Acti, equal to 0.17Acri at zero 
magnetic field and monotonically decreasing at stronger 
fields B. To avoid cumbersome expressions, we omit the 
terms with i = 3, 4 in our further analysis. 

The contribution to the Cooperon self energy due to 
the scattering off magnetic impurities is shown in Fig. 3. 
Details of our calculations arc presented in Appendix 
B. We find that the WL correction to the conductivity 
may be represented in the form of Eq. (30) with function 
Y(B,T) replaced by 

Y S (B,T)= [ %— t 1 (40) 

J 2Tcosh 2 e/2T v /(r(e)+ 7orb )T s ' 

— oo 

where the phase relaxation rate due to scattering off mag- 



netic impurities is 

r(£) = i_ (S 2 z ) + (S z )tznh(e + u s )/2T ~\ 1^ 
S(S + 1) T s 

provided that the impurity spin relaxation time is large, 
T T r(e - T) > 1. 

Spin correlation functions (S z ) and (§%) in Eq. (41) 
are defined by Eqs. (15). We notice that in strong fields 
us > T the impurity spins are polarized and the phase 
relaxation rate of electrons close to the Fermi surface 
|e| < T vanishes, in agreement with Ref. [8]. 

We consider the special case of magnetic impurities 
with S = 1/2 in more details. First, we discuss the 
meaning of the energy-dependent relaxation rate T(e), 
Eq. (41). In the S — 1/2 case, T(e) can be rewritten in 
the form: 

4 

T(e) = — \pi{l - n(s + u s )) + pM s + us)} , (42) 
3t s 

where Pf(l) = (2 coshcj5/2T)~ 1 exp (±us/2T) is the 
probability for the spin impurity to be parallel (an- 
tiparallcl) to the direction of the magnetic field and 
n(e) = [1 + exp(e/T)] _1 is the Fermi occupation num- 
ber for electrons with energy e at temperature T. We 
interpret Eq. (42) in the following way. Two processes 
contribute to the electron phase relaxation rate: i) an 
electron with spin up and energy e is scattered by an 
impurity with spin down to the electron state with spin 
down and energy e + us', ii) a hole with spin up and 
energy e is scattered by an impurity with spin up to 
the electron state with spin down and energy e + us- 
The probabilities of these processes are determined by 
the first and second terms in Eq. (42), respectively. 

In the limit of short impurity spin relaxation time, 
T T T(e ~ T) > 1, the function Yy 2 (B, T) has the form 
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+ 00 

,(B,T) = J - 



dx 
cosh 2 x 



2 cosh x 

3 cosh uj s /2T cosh(a; + uj s /2T) 



-1/2 



H~ 'YorbT's 



(43) 



We illustrate the behavior of function Yi/ 2 (J3, T) in the 
inset of Fig. 2 at temperature T = 0.3T o , where T is 
defined by Eq. (35) with S = 1/2. 

If the orbital effect of the magnetic field is small, 7 orb = 
0, we have 



F 1/2 (B,T) = 



sinh(3ws/4T) 



3 tanh(w 5 /2T) v /cosh(w s /2T) ' 



(44) 



and the weak localization correction to the conductiv- 
ity is given by Eq. (30) with Yi/ 2 (B,T) in the form of 
Eq. (44). In the opposite limit, when the orbital effect 
of the magnetic field on WL dominates over the effect 
of spin scattering, we may expand function Yi/ 2 (B, T) in 
(fsTorb) -1 and obtain 



^(B) = J°- ( 1 - 1^- . f J . (45) 
7r^Y 7or bV 3T s 7 orb rsmhws/Ty 

We notice that the second term in Eq. (45) has the same structure as the phase relaxation rate suggested in Ref. [18] 
to describe the WL correction to the conductivity. Therefore, the expansion, presented in Eq. (45), establishes the 
conditions of applicability of the suggested formula. 18 

Expression for the phase relaxation rate in Eq. (41) was derived in the limit of long impurity spin relaxation time, 
ttT(s <~ T) ^> 1. In the opposite limit, the spins at different moments of time are not correlated, and we have to 
substitute (S z ) 2 for (S 2 ) in Eq. (41). Then, the function Yi/ 2 (B,T) acquires the form: 

Yy 2 (B, T) = J (l - I tanh 2 % - § tanh ( x+ %) tanh % + 7orb . s ) (46) 



We observe from Eq. (41), that at strong magnetic field 
the phase relaxation rate exponentially vanishes. Conse- 
quently, higher order terms in the exchange constant J 
may become important. The fourth order diagrams for 
the Cooperon phase relaxation rate are shown in Fig. 4 
and the corresponding analytical calculations, presented 
in Appendix C, yield 



T(e) = 



ir 2 T 2 



2ir n s T s u}g(S + l) 2 r s ' 



(47) 



We emphasize that Eq. (47) represents the J 4 contribu- 
tion to the electron phase relaxation rate, while Eq. (41) 
represents the J 2 contribution. The corresponding small 
parameter of the expansion in powers of J 2 may be writ- 
ten as u/(n s T s ). 

As the applied magnetic field increases, both the J 2 
and J 4 contributions may become comparable, since the 
J 2 contribution decreases exponentially, while the J 4 
contribution at energy |e| ~ T decreases much slower, 
cf. Eqs. (41) and (47). The small factor v/(n s T s ) is com- 
pensated at strong enough magnetic field which we 



estimate from exp(— oos/T) = v/(n s T s ). We obtain 



gUB n s T s 



(48) 



At small fields, B < the relaxation rate decreases 
exponentially with the increase of B, whereas at higher 
fields this dependence is replaced by a slower power law 
r(e ~ T) oc T 2 /uj 2 . The counterpart of Eq. (45) at 
B > B* reads 



1 e 2 

Aai{ B) = --- 




3 n s r s T s 7 orb (5 + l) 2 w| 



• (49) 



V. CONDUCTANCE FLUCTUATIONS 

In this Section we study the effect of the impurity spin 
polarization on conductance fluctuations of a metal ring, 
weakly connected to the leads. 28 Magnetic flux $ pierc- 
ing the ring changes wave functions of electrons in the 
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T I 



FIG. 4: The Cooperon self energy diagrams to the fourth 
order in the exchange constant J. The spin-spin correlation 
function (wavy line) acquires correction due to the electron 
scattering off the spin. 



ring and, consequently, the ring conductance. The fluc- 
tuations of the conductance are usually characterized by 
the correlation function /C(At, A<&): 



£(At,A$) = «<7*(0)<7$+A$(Ai)})$, 



(50) 



where $ = irR 2 H is the magnetic flux through the ring of 
radius R, and At is the time lapse between the measure- 
ments of g$ and fM>+A$- The d.c. conductance, Eq. (50), 
is defined in terms of the current I(t), averaged over mea- 
surement time T m 



9<s>{t) 



rt+T m 
Jt 



01* (t) 



dV 



(51) 



v=o 



where V is the applied bias. During the measurement 
time T m the magnetic flux should be constant, otherwise 
the measured conductance is already averaged over dif- 
ferent realizations. 

We assume that the system satisfies the er- 
godic hypothesis and the averaging over the mag- 
netic field is equivalent to averaging over impurity 
configurations. 2,29 To calculate the conductance corre- 
lation function, Eq. (50), we apply the conventional av- 
eraging technique. 26 We consider only the case of short 
relaxation time tt <C At. In this case the spin correla- 
tion functions in Eqs. (16) should be taken in the limit 
t — > oo, so that x_l(t — > oo) = and Xz( T — * °°) = 
(S,) 2 . 

The conductance correlation function contains all har- 



monics 



£(A$) = ^ £ fc cos ( 2nk— I 

k=0 



*0 / 



(52) 



The amplitudes Kk consist of two parts, originating 
from the fluctuations of the diffusion coefficient, Fig. 5a, 
and the fluctuations of the electron density of states, 
Fig. 5b: 



_ e 4 D 2 f dede' , _ 2 e , -2 £ ' 

K,u = a- -tt — 7/ cosh — cosh — 

h (2irh) 2 i? 4 J 16T 2 2T 2T 



x ]T (4* W) + 4* W)) , 



(53) 



see Appendix D. Here a is a geometry dependent dimen- 
sionless factor of order of unity; functions A^\s,s') and 
Bf!\e,e') are defined by 



A ( ^(e,e')=2j T>i(e,e', -|) cos(2irkx)dx, (54a) 
fiW ( e , e ')=J Re \v 2 (e, e', ^) } cos{2Trkx)dx (54b) 



in terms of the diffuson components: 
T>i(e,e',q) = - 



i{e-e'+^e z ) + Dq 2 + lv + T l {e 1 e'Y 



(55) 



The index i in Eqs. (54) runs over different spin config- 
urations of the diffuson T>i(e,e' ,q), which are related to 
the common classification of spin wave functions of two 
spin- 1/2 particles in terms of triplet and singlet states, 
see Table II. In our notations, modes i = 1, 2 correspond 
to total spin = 1 with non-zero projections on the 
magnetic field = ±1, mode i = 3 represents total 
spin Sd = 1 with zero projection on the magnetic field, 
mj = 0, and mode i = 4 is a singlet spin configuration. 
Energies q£z represent the effect of the Zeeman splitting 
of conduction electron states on various diffuson modes. 
Coefficients q are given in Table II. 

We explicitly separated two additive components, 
Ti(e,e') and 7^, to the diffuson decay rate in Eq. (55). 
The component Ti(e,e') corresponds to the contribu- 
tions from spin-orbit interaction and scattering off mag- 
netic impurities. The component j v takes into account 
other processes, such as electron escape through the leads 
and scattering caused by electron-phonon and electron- 
electron interactions. 22 



A. Effect of classical spin impurities 

Following the spirit of Section IV, we first analyze con- 
ductance fluctuations in the presence of classical S S> 1 
spins. This analysis allows us to explore the effect of the 
impurity spin polarization at various relations between 
parameters of the system. In the next subsection we cal- 
culate the amplitude of conductance correlation function, 
Eqs. (54), for spins S <~ 1 in a metal with strong spin or- 
bit interaction. 

We consider systems with fast spin relaxation time 
tt -C At, where At is the time lapse between the current 
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* 0=00] 

FIG. 5: Diagrams a) and b) give the main contribution to the conductance correlation function. Diagram a) originates from 
the fluctuations of the diffusion coefficient and decreases only as 1/T at high temperature. Diagram b) is referred to as the 
contribution from the fluctuations of the electron density of states and exponentially vanishes at high temperature, see 28 . The 
Hikami box (black square) is the sum of two diagrams, shown in Fig. c). 



TABLE II: The decay rates of various components of the dif- 
fuson in a metal without spin orbit interaction. The compo- 
nents with i = 1,2, 3 represent a spin states with total spin 
one, Sd = 1. These spin configurations form a triplet in an 
isotropic system. Anisotropy due to the applied magnetic field 
lifts the degeneracy. The i = 4 component is a singlet spin 
configuration, Sd = 0. The last column shows coefficients q 
for the Zeeman splitting of conduction electrons. 



i 


|<Sd,md> 


Ti{H)T. 


?i 


1 

2 
3 
4 


Sd = 1, m d — +1 
Sd = l,m d = -1 
Sd = 1, m d = 
Sd = 0, md = 


i + (s.y/sr 

1 + (S Z ) 2 /S 2 
1 - (S z f/S 2 

i - <s z ) 2 /s 2 


+i 
-i 






measurements, see Eq. (50). In this case the solution 
of the diffuson equation, Eq. (A4), is described by the 
following energy-independent decay rates: 

Ti, 2 = (1 + (S z ) 2 /S 2 ) - + Uso, (56a) 

r 3 = (1 - (S z ) 2 /S 2 ) - + | 7 so, (56b) 

T 4 - (l-(S z ) 2 /S 2 )-. (56c) 

Here (S z ) is defined by Eq. (29a) and depends only on 
the ratio of the magnetic field and temperature B/T. 

The i = 1, 2 modes represent the interference of elec- 
tron states with opposite spin orientations. The exchange 
field of magnetic impurities produces different (opposite) 
contributions to the phases of these two states. Because 
the phase contributions for different electron trajectories 
fluctuate, the interference of electron states with opposite 
spins is suppressed even if all spins of magnetic impuri- 
ties are fully polarized. Thus, the polarization of impu- 
rity spins does not suppress the effect of electron phase 
relaxation due to the scattering off magnetic impurities 
for diffuson modes i = 1,2. According to Eq. (56a), the 



diffuson decay rates I\ 2 actually increase as the applied 
magnetic field increases. 

On the other hand, the i = 3,4 modes stand for the 
interference of two electron states with parallel spins. At 
strong magnetic field B 3> T/gfi^S, when all spins are 
polarized, the scattering off magnetic impurities provides 
equal phase shifts to both states and does not affect dif- 
fuson relaxation rates r 3:4 . That is why the contribution 
t° ^3Ai caused by the scattering off magnetic impuri- 
ties, vanishes as impurity spins become polarized, see 
Eqs. (56b) and (56c). 

Substituting the diffuson decay times from Eqs. (56) 
into Eqs. (53) - (55), we can describe the harmonics of 
the conductance correlation function Eq. (52) at an ar- 
bitrary value of magnetic field, ranging from B = to 
B 3> T/g/iftS. To analyze this crossover, we first con- 
sider a metal without spin-orbit interaction, 7 so = 0, and 
evaluate the integrals over energies e and e' in Eq. (53) 
for two limiting cases of low (T <C I\ + j v ) or high tem- 
perature (T >• Ti + 7 V ). 

In the low temperature limit we have: 28 

with 

Li = Jj^— (58) 

For simplicity, we omitted the Zeeman splitting of the 
conduction electron states, which actually modifies the 
i = 1,2 terms in Eq. (57). We emphasize that at low 
temperature both diagrams in Fig. 5a and 5b contribute 
to the conductance correlation function. 

In the limit of high temperature T 3> + j v , 
the contribution of the diagram in Fig. 5b is small as 
exp[-(2 7 r) 3 / 2 fc J R(T/L>) 1 / 2 ] and may be disregarded. The 
contribution due to fluctuations of the diffusion coeffi- 
cient, diagram in Fig. 5a, decays only as 1/T at high 
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temperature T: 



^- a 3 (27r?i)2^ / v2T 



<«£z\ L T Li 



-2-KkR/Li 



(59) 



where = \/D/T is the thermal length, and coeffi- 
cients q are presented in Table II. The function 



/(*) = 3- 



: cosh z — sinh z 
sinh 3 z 



(60) 



takes into account the Zeeman splitting; /(0) = 1 and 
/(z) w 12(z- l)e- 2z for z > 1. 

Regardless of scattering off magnetic impurities, the 
Zeeman splitting destroys the contribution of modes i = 
1, 2 to the conductance correlation function JC at e z ^> T, 
see Eq. (59). From Eq. (59) we conclude that at strong 
magnetic field, £z > T and Sus S> T, the amplitude of 
conductance fluctuations is no longer determined by the 
scattering rate off magnetic impurities, l/r s . Since modes 
i = 1,2 are also suppressed at ez > T even in a metal 
without spin impurities, we conclude that a strong mag- 
netic field restores the amplitude of conductance fluctua- 
tions up to the amplitude in the same ring as if magnetic 
impurities were absent. 

We also discuss another effect of the Zeeman splitting 
on the conductance correlation function. If the exchange 
constant J is antiferromagnetic, J > 0, the Zeeman split- 
ting vanishes not only at B — but also at some finite B, 
see Eq. (19), and thus produces another wrinkle in the 
dependence of ICk on the applied magnetic field. This 
non-monotonic behavior of the amplitude of the conduc- 
tance oscillations is reminiscent to the reentrance effect 
in superconductors. 23 

To illustrate the effects of the Zeeman splitting ez, 
we plot amplitude JC\ of the principal harmonic, fc = 1, 
as a function of B in Fig. 6. We choose the following 
values of the system parameters: 7 y r s = 2 and R = 
y/D/'Yp. The dashed line corresponds to the case of equal 
g-factors for the conduction electrons and impurity spins, 
i.e., ez = u)s, and negligibly small exchange contribution 
in Eq. (19), n s J(S) < u s . The limit of e z > T, when 
only terms i = 3,4 in Eq. (59) survive, is represented by 
the solid line. Finally, the dotted line demonstrates the 
reentrance effect due to the antiferromagnetic impurities 
at some specific value of n s J = 3.5T with T being the 
temperature. 

Now we notice, see Eqs. (56), that the spin-orbit in- 
teraction suppresses the contribution to the amplitude of 
the conductance fluctuations, originating from the diffu- 
son modes with i = 1, 2, 3 (5<j = 1). In the limit of strong 
spin-orbit interaction (7 so » D/R 2 ) only the contribu- 
tion from the singlet (Sd — and i = 4) mode survives: 



JCk = a 



3 (27rft) 2 R 3 



L T La c -2irkR/L 4 



(61) 




with length L4 defined by Eqs. (56c) and (58). The 
function ICk for a metal with strong spin orbit scattering 



4 6 

g^SB/T 

FIG. 6: Amplitude JCi(B) as a function of the applied mag- 
netic field B for a metal with classical (S ^> 1) magnetic 
impurities in the absence of spin orbit interaction. A non- 
monotonic behavior of ICi due to the Zeeman splitting of 
the conduction electron states is illustrated for two cases: i) 
ez = (dashed line), ii) ez = us — 7T(S Z ) (dotted line), 
here the second term represents the average exchange field 
2n s J{S z ) with n s J = 3.5T. The limit ez —> 00 is presented 
by a solid line. We choose 7 v t s = 2 and R = \J Dj^ v . 



monotonically increases as a function of B/T, provided 
-f v remains constant. At strong magnetic field, Sujs S> T, 
magnetic impurities do not affect the conductance fluc- 
tuations. 

We emphasize that the non-monotonic behavior of the 
amplitudes JCk of the conductance correlation function 
originates from the diffuson modes with i = 1,2, which 
do not contribute to the correlation function in metals 
with strong spin-orbit interaction. Thus, the monotonic 
or non-monotonic behavior of the amplitudes ICk distin- 
guishes materials with or without spin orbit coupling. 



B. Effect of quantum spin impurities 

In this section we perform a quantum calculation to an- 
alyze the effect of polarization of impurity spins with S ~ 
1, when the semiclassical description of electron scatter- 
ing off magnetic impurities is not applicable. We consider 
metals with strong spin orbit interaction, 7 so 3> D/R 2 , 
so that only the singlet component of the diffuson sur- 
vives. In this case the calculation of the diffuson self en- 
ergy is similar to the calculation of the proper quantity 
for the Cooperon, see Sec. IV C. Details of calculations 
and the expression for the amplitude of the conductance 
fluctuations are presented in Appendix E. 

Below we focus on the high temperature case, T 3> 
r(e,e' ~ T), which is reached at sufficiently strong mag- 
netic field even if t s T « 1. At high temperature, the 
integral over the difference of e and e' converges fast, 
\e — e'\ < T(e, e' ~ T). This observation allows us to per- 
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FIG. 7: The dependence of conductance oscillations JCi(B) as 
a function of the applied magnetic field B for four values of 
impurity spins: S = 1/2 (solid line), S — 3/2 (dashed line), 
S = 5/2 (dotted line), S = 7/2 (dash-dotted line). We choose 
■y v T s = 1/2 and R = sjD/^^. 



form integration over the energy difference. We obtain: 



a- 



DL 2 T 



8 {2nh) 2 i? 3 

27vkR y /r(e)+ rf /VD d£ / T 



I 



y/T(e) + 7v cosh 4 e/2T' 
with the dcphasing rate in the form of 



T(s) = 



1 - 



(^) 2 + (^)tanh(g + u s )/2r 
S(S + 1) 



(62) 



(63) 



We represent the corresponding curves of ICi as a function 
of the applied magnetic field B in Fig. 7 for several values 
of S, assuming t s T 3> 1 and for the following values of 
the system parameters: 7 v t s = 1/2 and R = yj Djj^. 
The shapes of the curves corresponding to various val- 
ues of the impurity spin S are different from each other. 
Particularly, we conclude that the conductance fluctua- 
tions are faster restored by magnetic field in metals with 
a larger value of S. The fitting of experimental results 
by the curves JCk(B) given by Eq. (62), may provide the 
impurity spin parameters, such as its value S and the 
gyromagnetic factor g. 

As the polarization (S z ) and the Zeeman splitting us 
increase, the diffuson decay rate r(e ~ T) to the second 
order in the exchange constant J becomes exponentially 
small, see Eq. (63). In this case we have to take into 
account higher-order contribution in J, which is repre- 
sented by Eq. (47). The crossover from the phase relax- 
ation rate in the form of Eq. (63) to the form of Eq. (47) 
takes place at magnetic field B <~ £?*(T), see Eq. (8). 



VI. CONCLUSIONS 

The main question addressed in this paper, is how the 
polarization of magnetic impurities in a metal affects 



the interference contribution to its conductance. It is 
well-known that at weak magnetic field, which does not 
cause such polarization, the electron scattering off local- 
ized spins results in suppression of the weak localization 
correction to the conductivity and in suppression of the 
mesoscopic conductance fluctuations. If these interfer- 
ence contributions can be restored by an application of a 
stronger spin-polarizing field, then a transport measure- 
ment may serve as a test for the presence of magnetic 
impurities in a sample. This possibility was the main 
motivation of the presented study. 

We obtained analytical results for the WL correction 
to the conductivity and for the amplitude of conductance 
fluctuations valid for an arbitrary magnetic field. We 
found that the conditions for the restoration of weak lo- 
calization correction are quite stringent. The weak lo- 
calization correction can be substantially enhanced by 
an application of a magnetic field only in samples of very 
small size (to avoid the orbital effect of the magnetic field) 
and made of a light-element material, to make the spin- 
orbit scattering negligible. If these conditions are met, 
the spin polarization may completely eliminate the effect 
of magnetic impurities. At intermediate fields B, the 
weak localization correction depends on the ratio B/T, 
and we find the corresponding crossover function span- 
ning the full range of spin polarizations. 

The effect of spin polarization on the mesoscopic con- 
ductance fluctuations is much more robust. The ampli- 
tude of fluctuations is restored if variations of a spin con- 
figuration in time arc suppressed by the external mag- 
netic field. At the same time, there is no detrimental 
orbital influence of the applied magnetic field on meso- 
scopic fluctuations. The spin-orbit interaction also docs 
not spoil the effect of spin polarization on the amplitude 
of fluctuations. We evaluated the amplitude of the "ft/e" 
Aharonov-Bohm oscillations of the conductance correla- 
tion function. In a metal with strong spin orbit interac- 
tion, the polarization of magnetic impurities restores this 
amplitude up to its value characteristic for the host ma- 
terial in the absence of magnetic impurities. We find the 
full crossover function, which describes the amplitude of 
fluctuations at arbitrary value of B/T. The use of our 
results for the analysis of transport measurements may, 
in principle, yield such characteristics of magnetic impu- 
rities as the values of their spin and g-i actor, which are 
hard to access directly. 
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APPENDIX A: COOPERON AND DIFFUSON IN 
A METAL WITH CLASSICAL MAGNETIC 
IMPURITIES 



ments, we use the Cooperon which is formally defined as 
an average of electron Green's functions 



In calculations of the weak localization correction to 
the conductivity of a metal with classical magnetic mo- 



J 



C 



^(tJ_*L,^_^ >t + +tr;P y)j (t + +tr _4_^- ) = {G W( t + 4;r ) r')G^(i r ./7:v.r'.) (Al) 



After the standard procedure 1 we obtain the following expression for the Cooperon: 

d 



^ + ^q 2 ) i + n c Q +n c s (r) 



C(r,r',r;q) = l(5(r-r'), 



(A2a) 



where Hq is a matrix in spin space representing the Zccman splitting of conduction electron states and spin-orbit 
interaction: 
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and H^ ms describes scattering off magnetic impurities: 



/(1-X**(r)) 



K(r) = - 
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(A2c) 



The solution of Eqs. A2 yields the Cooperon in the form of Eq. (21). 

The diffuson appears in this paper in calculations of the conductance correlation function and is defined as 



V c 



j6 \ 2 7 2 
It satisfies the following equation: 



,t+-t-;r,r')S(t+-t--t++t-) = (G™ (t+, t+; r, r')G^ } (tj.tj"; r', r)). (A3) 



( W + J Dq 2 )l+^+Wf(T) 2? q (w,T) = l, 



(A4a) 



where the matrix 7t^ docs not depend on time variable r and describes the Zeeman splitting of conduction electron 
states and spin-orbit interaction: 
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and n^(r) represents scattering off magnetic impurities: 



















2 












2 


!) 















(A4b) 



/(1-X«W) 



V -xl(r) 





(i + x«W) o 

+ 











-X-l(t) 



(1-X«W) 



(A4c) 



r 



Solving Eqs. (A4), we obtain the diffuson in the form APPENDIX B: COOPERON DECAY RATE DUE 
of Eq. (55). TO SCATTERING OFF QUANTUM MAGNETIC 

IMPURITIES 



In this Appendix we derive the Cooperon decay rate 
in a metal with quantum (S <~ 1) impurity spins. The 
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FIG. 8: The Keldysh contour allows to take into account 
time ordering of quantum operators. The time subscripts if,b 
mean that the operator is taken at time t at the forward (tf) 
or backward (t b ) part of the Keldysh contour. 



Hamiltonian of interaction of conduction electrons with 
magnetic impurities, Eq. (12), has the form: 

H m = J + S-&+ + (Bl) 

where a x , a y and a z are the Pauli matrices and o± = 
(<j x ± iby) /2. We notice that the first term in Eq. (Bl) 
represent scattering without spin flip, while two other 
terms correspond to spin flip scattering. 

We will follow the standard Keldysh formalism: first, 
we define electron and spin operators on the forward and 
backward parts of the Keldysh contour, see Fig. 8, and 
then perform rotation in the Keldysh space. As the re- 
sult, we obtain the standard form of the electron Green's 
function in terms of the retarded G^(e,p), advanced 
G*( A ) (e, p) and Keldysh (e, p) components. The spin 
operator in the rotated basis acquires the form 

s * {t) ~ 2 ^ sm sm ) ' (B2) 

where 5, s ' a (i) — Si(tf)± Si(tb)- In our calculations of the 
impurity average conductivity, we work in the Born ap- 
proximation, therefore the scattering is completely char- 



acterized by the bilinear spin correlators: 

= \(SmS s 3 (t)), (B3a) 
- \{St(t)Sf{0)), (B3b) 

c-m = \{smsm)' ( B3c ) 

which are counterparts of the Keldysh Green's function 
of an interaction field, see 30 . 

For the Zeeman Hamiltonian H sp i D — g[iBBS z , spin 
correlator components with i = j = z have only the 
Keldysh clement L\ z {t): 

& zz {t) = {Sl{t)S s M) = <5f), (B4) 

and £ T zz (t) vanish due to the commutation relations of 
the spin operator along z axis and the Hamiltonian H sp i D . 

Now we present correlators of the spin components per- 
pendicular to the magnetic field: 

Ct+(t) = Oi-Qe**"* (S x ) , (B5a) 
CL + {t) = -Stfe*"'* (S t ) , (B5b) 
C k _ + {t) = e*"'* (S(S + 1) - (5f)) . (B5c) 

Using the formal definitions of Eqs. (B4) and (B5), we 
calculate the Cooperon decay rate, which has the mean- 
ing of the Cooperon self energy, see Eq. (21). The dia- 
grams which contribute to the decay rates of Cooperon 
modes i = 1,2 are shown in Fig. 3a and 3b and corre- 
spond to the electron self-energy due to magnetic impu- 
rities: 



Er(£) = /(5p/ ^{(^H+^H) Gr ( e + ^P) + ( £ - + H + ^M)G k (£ + c,,p)}, (B6a) 
^ = J jSj^J ^{(^H+^H) Ga ( £ + ^P) + ( £a - + H+^H)G k (e + ^,p)}. (B6b) 

I 

The elastic part of the self-energy is expressed only in r a b(e) = m s J 2 [S a (e) — S b (e)] and can be written as 
terms of the retarded and advanced Green's function 
and is not affected by the electron distribution, while 1 
the inelastic part of the self energy contains also the 

r ab ( £ ) = - 

Keldysh component of the electron Green's function. 

Consequently, for inelastic processes the electron distri- The firgt term in Eq (B?) originates from terms in 

bution is important, since such processes should satisfy Eqg (B6b)j containing the Keldysh components of the 

the fermion exclusion principle. The contribution to the gpin correlators r k (w) and fr_Au}). We used the prop- 

Cooperon decay rate from diagrams in Figs. 5a and 5b is erty of ^ corrdators £ k^ = Q) + + = Q) = 

5(5+1), arising from their definitions, see Eqs. (B3). We 
observe that the first term in Eq. (B7) is the total electron 



tanh 



: + ui Im {£i + (w)} dw 
~2T 5(5+1) 2T 



(B7) 
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scattering rate on magnetic impurities l/r s . The second 
term represents the rate of inelastic scattering processes, 
accompanied by the energy transfer. Particularly this 
term explicitly contains the electron distribution func- 
tion in the form 1 — 2n(e) = tanhe/2T, which takes into 
account the Pauli exclusion principle for electrons. 

We identify the diagram in Fig. 3c with the vertex cor- 
rection to the Cooperon self energy. For the Cooperon 
component with parallel electron spins, only the com- 
ponent of the scattering on magnetic impurities without 
spin flip remains: 



r c (e) = 



(S 2 Z ) 



T s 5(5+1) 



(B8) 



The decay rate of Cooperon modes i = 1,2 is given by 
the sum of r a b(e) and r c (e): 



m = - 



(S 2 ) + (S z ) tanh(e + to s )/2T 



5(5 + 1) 



• (B9) 



When deriving Eq. (B9) we assumed that the conduc- 
tion electron spins and average impurity spins are paral- 
lel, so that uj s > and (S z ) > 0. The opposite case is 
also described by Eq. (B9) with uj s < and (S z ) < 0. 
Equation (B9) determines the energy-dependent phase 
relaxation rate in the WL correction to conductivity, see 
Eq. (30). 



APPENDIX C: THE /-ORDER CONTRIBUTION 
TO THE ELECTRON SELF ENERGY 

In this Appendix we derive the correction to the 
Cooperon self energy to the fourth order in the exchange 
constant J. For this purpose we will use the rules of 
the Keldysh diagrammatic technique, outlined in Section 
IV C. 

The contribution to the Cooperon self energy is shown 
in Fig. 4. The electron-hole loops can be treated as 
the self energy g of the spin-spin Keldysh Green's func- 
tion, responsible for the spin relaxation as a result of 
the electron inelastic scattering. The imaginary part of 
this electron-hole loop determines the Korringa spin re- 
laxation rate. 7 In the Keldysh formalism, we can distin- 
guish three components of the self energy, related to the 
advanced, retarded and Keldysh elements of the spin-spin 
Green's function. 

We have for the difference of retarded and advanced 
components of the spin-spin self energy the following ex- 
pression 



2-kv 2 J 2 uj 1 



and the Keldysh component is 



g k H = coth^(^H- e »). 



(Cl) 



(C2) 



This result is of no surprise. Indeed, the self energy is 
not related to the type of interaction, since it is uniquely 
determined by the electron-hole loop, rather than the 
interaction lines. 

Here we are interested in the strong magnetic field 
limit, us S> T, so that only two impurity spin states with 
S z = 5 and S z = 5 — 1 arc involved. Under this condi- 
tion, the diagrams in Fig. 4 can be represented in terms 
of the modified spin-spin correlators, which to the second 
order in the exchange constant v 2 J 2 have the form: 

C\\u) = 4°- r '>)e r '>)4°-' a M ( C3 ) 



#M = 4°1>)2V)£^>) 



(0)a 



(C4) 



where C^ T ' a {u3) = (S)/(u> - ± iO) are the bare spin- 
spin correlation functions, defined by Eq. (B5). Using 
this notations, we reduce diagrams in Fig. 4 to diagrams 
in Figs. 3a and 3b. Then, after the standard calculations, 
we obtain 



2 T 2 , 2 

r 4 (e) - 2^ 3 J 4 5 2 — 7 . 



(C5) 



Using the definition of t s , we rewrite Eq. (C5) in the form 
of Eq. (47). 



APPENDIX D: CONDUCTANCE CORRELATION 
FUNCTION 

In this Appendix we outline the calculation of the con- 
ductance correlation function /C(A$), see Eq. (52). The 
contributions to /C(A$) originate from two diagrams, 
shown in Fig. 5. We have 

£(A$) = /C (a) (A$) +/C (b) (A$), (Dl) 

where the first term /C^(A$) is the contribution from 
fluctuations of the diffusion coefficient 



/C (a) (A$) = a- 



3 4 D 2 f dede' 



1 



(2irh) 2 R 4 J 16T 2 cosh 2 e/2TcoshV/2T 



+ 00 



^ EE 



V, 



A$ 



(D2) 



and second term in Eq. (Dl) originates from fluctuations 
of the electron density of states: 



dede' 



16T 2 C osh 2 e/2T coshV/2T 
A$ 



(D3) 



Here T>i(e,e',q) is the ith diffuson component and the 
sum over n runs through all discreet values of momentum 
q n = n/R. We apply the Poisson summation formula 



E /(»)= E //(»)< 



2-nikn 



dn 
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to Eqs. (D2) and (D3) and obtain Eq. (52) with ampli- 
tudes JCk denned by Eqs. (53) and (54). 



we obtain the following expression for the self energy of 
electron Green functions, see Eqs. (B6b): 



APPENDIX E: DIFFUSON DECAY RATE DUE 
TO SCATTERING OFF QUANTUM MAGNETIC 
IMPURITIES 

We consider the singlet contribution to the conduc- 
tance fluctuations in the limit of short spin relaxation 
time 7T, so that the components of spin perpendicu- 
lar to magnetic field arc not correlated at the char- 
acteristic time At between conductance measurements, 
f±(At) — > in Eqs. (16). In this case the diffuson self- 
energy is defined by diagrams similar to ones shown in 
Fig. (3), but with the opposite direction of the advanced 
Green's function. 

The diffuson self energy consists of three terms; two 
of them coincide with the self energies S r (e) and S a (e) 
of retarded and advanced electron Green functions, re- 
spectively, and the third term is the vertex correction 
S v (e). Using the formalism, presented in Appendix B, 



£ r ' a (e) = w- 

ZT S 



1 (5 z ) 2 + (5 z )tanh( £ + ^)/2T 



5(5 + 1) 



± i 



. (S z ) f t&nh{e + uj)/2T du 
'5(5+1) J u + u s 2^ 



• (El) 



The vertex correction to the diffuson originates from the 
impurity spin correlator, taken at time difference At, see 
Eq. (50). Since we assume that the impurity spin relax- 
ation time is much smaller than the delay time At, the 
spin correlator is \ z — (5 Z ) 2 and we have 



^{e) = -{S z )\ 



(E2) 



Combining the electron self energy part and the vertex 
correction given by Eqs. (El) and (E2), we obtain the 
diffuson self energy in the form: 



r(e,e') = 



1 - 



+ i 



<*>/ 



(5 Z ) 2 (S z ) ( _ 

5(5+1) 25(5+1) V 2T 

tanh(e + u)/2T - tanh(e' + uj)/2T duj 



w s e' +uj, 

h tanh — — 

2T 



LU + LU S 



2tt 



(E3) 



For the calculation of the conductance fluctuations, we have to consider a finite difference between the energies e and 
e', transferred along the advanced and retarded Green's functions. The last term in Eq. (E3) originates from the real 
part of the electron self energy and vanishes if e = e' . This term represent renormalization of the electron density of 
states. 

The amplitude of the fcth harmonic of the conductance correlation function is expressed in terms of and , 
see Eqs. (53) and (54): 



4 4) 



4 4) = 



2tt 



R 2 e 



8^2 



Rc- 



-Imr(e,e') 
1 + 2nk 



Im 



exp (-27rfciV7 v + r ( e > e + *( e - e')/VDj 
^ lip + T(e,e')+i(e-e') 

RVD \ ( „ , ^l v +T{e,e')+i{e-e' 



^ v +T(e,s')+i(e-e') / 
I 



exp —2irkR- 



(E4a) 
(E4b) 



At high temperature Tt s » 1 integral in Eq. (53) over and (63). 
e — e' converges at \e — s'\ < l/r s , and we obtain Eqs. (62) 
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